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NOMENCLATURE 


c 


m 


K/c 


P 

P 

oo 

Hoc 

T 

V 


x/c 

z/c 


Airfoil chord length, m 
Drag coefficient, D/qooC 
Lift coefficient, L/q^c 

Moment coefficient about the quarter chord, 
M/q c 2 

Pressure coefficient, (P - P 00 )/q 00 

Roughness height 

Local static pressure, N/m 2 

Free stream static pressure, N/m 2 

Free stream dynamic pressure, N/m 2 

Temperature, °F 

Velocity in knots 

Horizontal coordinate 

Vertical coordinate 


a, AOA Angle of attack, degrees 
Flap deflection, degrees 


in 




INTRODUCTION 


The test program described in this report is an extension of 
a study begun in 1981 to provide needed information on the perfor- 
mance degradation of airfoil sections resulting from rime and glaze 
ice accretions . Its primary objectives were: 

1) To expand the current database of performance data on the 
63A415 with simulated ice to include flap deflection. 

2) To further study the flowfield in the area of the ice 
accretion through pressure distributions and flow visualization 
techniques, which can then be used to evaluate the accuracy of the 
theoretical analysis methods currently being developed. 

3) To obtain data on a simulated glaze ice shape that scales 
down to a 6 inch chord model and will be tested in the OSU Tran- 
sonic Airfoil Wind Tunnel Facility. These data will be used to 
compare the aerodynamic qualities of the NASA Icing Research Tunnel 
and the OSU tunnel, and to evaluate a lift measuring system based 
on wall pressures. 

Mr. Richard Freuler, Senior Computer Specialist at the 
Aeronautical and Astronautical Research Laboratory, developed the 
software needed for the data acquisition system. Mr. Steven 
Thompson, an under-graduate research assistant, modified the 
software and performed the data reduction for this test. 
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EXPERIMENTAL METHOD 


Equipment 

Testing was performed in the NASA Lewis 6' x 9' Icing Research 
Tunnel (IRT) . The airfoil model used was the NACA 632~A415 with 
a 1.36 m chord and a moveable flap with deflections of 10°, 20°, 
and 30°. The airfoil and flap were pressure tapped using 1/8" OD 
strip- a- tube attached to the airfoil surface. In addition, the 
model was fitted with five simulated ice shapes (Figures 1A-1E) : 

1) Generic Glaze 

2) Glaze 3° 

3) Rime 3° 

4) Glaze 7° 

5) Rime 7° 

Aerodynamic data were taken on the first three shapes and flow 
visualization was performed on all five. The Generic Glaze shape 
was derived from the work of Ingelman-Sundberg 1 . This shape was 
chosen because it scales to a convenient size on the 6" chord model 
which will be tested in the OSU Transonic Airfoil Wind Tunnel 
Facility. 

The Glaze 3, Rime 3, Glaze 7, and Rime 7 shapes were chosen 
from a series of ice growths generated during an actual ice accre- 
tion study in the NASA Lewis Icing Research Tunnel^. They repre- ^ 
sent typical climb , high angle of attack and low velocity , and 
cruise, low angle of attack and high velocity conditions. 

In order to add the surface roughness characteristic of natural 
ice shapes, aluminum oxide grit with a K/C = .00058 was attached 
to the glaze shapes with a spray acrylic adhesive. A grit with a 
K/C = .0012 was added to the rime shapes. 

On-line data acquisition and reduction were accomplished using 

O 

the OSU Digital Data Acquisition and Reduction System (DDARS - 
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figure 2) . The central processing unit is the DEC LSI-11 micro- 
computer. Input and output is through a teletype terminal and mass 
data storage through a twin floppy disc drive system. Analog data 
signals from the transducers and wake probe slidewire systems are 
fed into an analog front end which conditions the signal and con- 
verts it into a digital format. 

Airfoil pressures were obtained through a Scanivalve trans- 
ducer arrangement, while drag data were measured using a wake probe 
with total and static ports. The voltages from these systems as 
well as those from tunnel total and tunnel static transducers were 
input to the analog box and then to the computer for on-line reduc- 
tion (figure 3). 

In order to visualize the flow in the leading edge region, 
a splitter plate^ was constructed which could be positioned between 
the upper and lower segments of the attached ice shape. (See 
figure 4) . Small drops of oil-based paint were then applied to the 
plate in the regions of interest and the tunnel then brought up 
to speed. Videotape was made of the movement of the drops and 
still photographs were taken after no further movement was observed. 

Data Reduction 

The DDARS system provides the test engineer quick-look pres- 
sure distributions as well as integrated values of C^, Cm, and C^. 
This permits maximum use of tunnel time. 

An interactive computer program was written for the final data 
reduction on the OSU Harris/6 computer system. The raw data files 
from the IRT test were transferred to the Harris from the LSI-11 
microcomputer. The program converts Scanivalve voltage from each 
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model tap into a pressure coefficient. The user is given a plot of 
the final Cp distribution for each element (main and flap) on a 
Tektronix CRT and can control any re-reduction required using the 
terminal cursors. The program then integrates the distribution to 
get lift and moment coefficients. 

The drag coefficient is calculated using the Jones Equation 5 . 
The wake is displayed on the graphics terminal and the user enters 
the integration limits using cursors. If the operator sees that 
the probe traverse was not large enough to capture the full wake, 
that run reduction can be bypassed. 
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RESULTS AND DISCUSSION 


Aerodynamic Measurements 

Data were taken on the following simulated ice accretions 
as well as the clean airfoil; 

1) Rime 3 Rough 

2) Glaze 3 Rough 

3) Generic Glaze Smooth 

4) Generic Glaze Rough 

In addition, for each configuration flap deflection was varied 
from 0-30 degrees. 

The glaze ice C p distributions show the characteristic adverse 

pressure gradient where the flow is forced to negotiate the large 

change in surface slope at the tip of the horns. These pressure 

spikes promote separation and tend to decrease Co and increase 

max 

the drag coefficient. The separated zone is clearly seen as a 
region of constant pressure in the C p distribution in the area 
behind the glaze ice horn. 


From the pressure distributions, it is observed that the flap 
was stalled for most of the runs. This separation is again char- 
acterized by a region of constant C p . A previous investigation by 
W. R. Krolak^ on the Beechcraft Sundowner, equipped with a NACA 
63A415 airfoil, shows this same trend in flight test data. 


From Table I and figure 5, it is clear that the penalties 
associated with ice show up in reductions in C£ max and a gtall . 
The G3 shape showed a reduction in C£ max over the clean case of 
0.2 - 0.4, and a reduction in a g£all as muc h as 4° for the 
6 f = 30° case. Similar reductions were seen for the generic and 
rime shapes . 

Due to the position of the wake probe, drag data could only 
be taken on 6£ = 0° cases. Cumulative plots of C £ vs. C(j show 
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TABLE I 


PERFORMANCE DEGRADATION WITH SIMULATED ICE 



CLEAN 

03 

GEN 

R3 


10 

20 

30 

10 

20 

30 

10 

20 

30 

10 

20 

30 

c *max 

1.8 

2.0 

2.2 

1.4 

1.7 

2.0 

1.2 

1.5 

1.7 

1.5 

1.75 

1.95 

^srall 

14.0 

12.5 

11.5 

10.5 

0.5 

7.5 

7.5 

- 

5.5 

10.5 

8.5 

6.5 

'*1.0 

-6.5 

-10.0 

-13.0 

-6.0 

- 

-12.5 

-6.0 

- 

-10.0 

-6.0 

- 

-11.0 


the increase in drag caused by the ice shapes. For example, at 
Cj? . = -4, a 20% increase in drag over the clean airfoil was observed 
when the R3 rough shape was attached, and a 30% increase for the 
G3 rough shape. Interestingly, the presence of roughness on the 
Generic Glaze shape was not found to be very crucial. This is due 
to the large laminar separation bubble in the region of the ice 
shape, which tends to be the prominent source of pressure drag. 

From the cumulative plots of C m vs C^, it is observed that 
at the lower lift coefficients the effect of the ice shape is 
almost negligible. However, at the higher Cjj,'s, for example at 
low speed with the flap deflected, more positive Cm's were observed 
with the simulated ice than for the clean airfoil. 

Flow Visualization 

Using the splitter plate arrangement, discussed previously, 
the flow about the simulated shapes was recorded. Of particular 
interest were the separated zones observed with the glaze shapes. 
These laminar separation zones were photographed and later the 
coordinates of the separated streamline were digitized from these 
records. Figures 6 and 7 are representative of the observed flow 
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patterns . Figure 6 clearly shows the Generic Glaze shape at 
a = 1.7° with its separated zone behind the horn. Figure 6 is of 
the same configuration but at a = 5.6°, and clearly shows the 
characteristic recirculation region. Figure 8 shows the G3 shape 
at a = 5.6°. 

The authors discovered during the analysis of the photos that 
the splitter plate extended too far into the flow ahead of the 
stagnation region between the glaze ice horns. The splitter plate 
boundary layer then separated due to the adverse gradient 
from the airfoil flowfield. This 3-D flowfield created vortices 
which were shed downstream and affected the flow patterns recorded. 
This is particularly evident in figure 7 where the streamlines 
converge due to the influence of these shed vortices. However, 
qualitatively the data provides some interesting clues to the 
shape and extent of the laminar separation bubble. 

Further investigation was performed at Ohio State using two 
different splitter plate configurations. A scaled-down version 
of the splitter plate utilized in the Lewis IRT and a smaller one 
with the leading edge reduced were tested on a GAW-1 airfoil with 
a simulated ice shape. Flow visualization techniques confirmed 
the authors hypotheses that vortices were shed downstream due 
to the severe pressure gradient induced by the ice shape on the 
splitter plate. It was observed that the reattachment point was 
shortened by as much as 3% under these test conditions as a 
result of the larger splitter plate. This value cannot however, 
be directly applied to the 63A415 airfoil, in the Lewis test. 

Rather, the reader should realize that qualitatively this shows 
that the observed reattachment point was moved forward due to 
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the presence of the splitter plate. In addition, it must be pointed 
out that this method of visualization does not actually display the 
position of the separated streamline. Rather a position above the 
zero velocity line in the separated zone between the recirculating 
flow is measured. 

Presentation of Data 

A tabulated run summary is included in the appendix of this 
report. It is organized by configuration: 1) clean, 2) rime 3 
rough, 3) glaze 3 rough, 4) generic glaze rough and 5) generic 
glaze smooth. Following these tables are the cumulative plots 
of C. vs a, C vs C 0 , and C 0 vs C,. Lastly, the pressure distri- 
butions are included and ordered in the same sequence as the run 
summary tables. 

Data reported with zero flap deflection was taken at 
approximately Re = 4.2 x 10 and M = 0.13. Due to the large 
loads on the model, data at all flap deflection angles greater 
than zero, were taken at approximately Re = 3.3 x 10^ and M = 0.10. 
No tunnel wall corrections have been made in the data. 
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SUMMARY AND CONCLUS IONS 


A typical general aviation airfoil, the NACA 632-A415, was 
outfitted with simulated ice accretions and tested in the NASA 
Icing Research Tunnel. Pressure distributions were obtained for 
a variety of flap deflections and angles of attack. As a result 
of this study, the following observations can be made; 

1) The airfoils with simulated ice shapes showed large 
increases in drag and heavy penalties in C£ max and a sta n. A 
shift in otLo was also observed. These reductions in performance 
would be of particular importance to the pilot in a landing con- 
figuration with the flap deployed and power reduced. 

2) Measured pressure distributions and flow visualization 
show the separated zone behind the horn of the glaze shapes and 

the severe adverse pressure gradients which lead to the separa- 
tion. 

3) Surface roughness for the Generic Glaze shape was not a 
crucial factor in the drag observed. Rather, the prominent effect 
was the large separated zone. 

Further investigation is necessary to document the flow character- 
istics reported. More detailed pressure distributions should be 
obtained, particularly in the region behind and between the glaze 
ice horns. Also, while flow visualization provides valuable 
insight into the flow in the separated zones, quantitative data 
must be gathered here before an analytical model can be developed. 


i 
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FIGURE 1A. R3 ICE SHAPE 
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FIGURE IB. R7 ICE SHAPE 
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FIGURE 1C. GLAZE 3 SIMULATED ICE ACCRETION 
AND PRESSURE TAP LOCATIONS 
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FIGURE ID. GLAZE 7 SIMULATED ICE ACCRETION 
AND PRESSURE TAP LOCATIONS 
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FIGURE IE . GENERIC GLAZE SIMULATED ICE ACCRETION 
AND PRESSURE TAP 10 CATIONS 
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FIGURE 2. OSU DIGITAL DATA ACQUISITION 
AND REDUCTION SYSTEM 
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FIGURE 3. OSU DATA ACQUISITION SYSTEM 
AS USED IN THE NASA LEWIS IRT 
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FIGURE 4. 63A415 WING WITH SPLITTER 
IN LEWIS ICING RESEARCH TUNNEL 
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FIGURE 5. CHANGE IN C 


WITH SIMULATED ICE SHAPES 


£ 

o max 

(6^ = 0 Cases from 1982 IRT Test 2 ) 
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FIGURE 7. SPLITTER PLATE PHOTOGRAPH OF UPPER 
SURFACE OF THE GENERIC CLAZE ROUGH ICE SHAPE 





FIGURE 8. SPLITTER PLATE PHOTOGRAPH OF UPPER 
SURFACE OF THE G3 ROUGH ICE SHAPE 
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Run Summary 


RUN * 

AOA 

FLAP 

DEF 

V < KT > 

T °( F ) 

PRESS . 
ALT . ( FT ) 

CL 

CD 

CM 

CLEAN 

48 

10*6 

0.0 

103.43 

77.0 

838*0 

1.3393 

0*0256 

- 0.095 

49 

11*6 

0.0 

101,92 

75.0 

832.3 

1*4410 

- 0 . 108 

60 

61 

”5 * 4 

0*0 

102.23 

50.0 

852.9 

- 0.2326 

0.0110 

- 0*044 

-2 * 4 

0*0 

101.39 

52.0 

844.3 

0.0710 

0.0105 

- 0.046 

62 

63 

- 0*4 

0*0 

100.18 

69.0 

839,2 

0.3126 

0.0108 

- 0.049 

1*6 

0 . 0 

102.84 

75.0 

869.4 

0,5431 

0.0112 

- 0 , 059 

64 

65 

3*6 

0*0 

101.05 

67.0 

857.2 

0.7220 

0*0128 

- 0.059 

5*6 

0 . 0 

102.25 

70.0 

870.5 

0.9020 

0.0144 

- 0*059 

66 

7*6 

0.0 

101.64 

72.0 

869*9 

1.1002 

0.0169 

- 0.075 

67 

68 

8*6 

0*0 

101*03 

74*0 

871.1 

1.1691 

0.0196 

- 0 , 077 

9*6 

0*0 

103.57 

75.0 

892.1 

1.2279 

0.0220 

- 0.075 

69 

70 
81 
82 

10*6 

0*0 

101*74 

71.0 

878,5 

1.3198 

0.0275 

- 0*092 

11*6 

0.0 

101.25 

73.0 

880,2 

1.4225 

0.0300 

- 0.090 

-2 • 4 

0 . 0 

102.59 

62.0 

1083.5 

0*0881 

0*0118 

- 0.049 

- 0*4 

0 ♦ 0 

101.79 

57.0 

1079,6 

0,3017 

0.0124 

- 0.051 

83 

1*6 

0.0 

102*24 

68*0 

1087,2 

0*5059 

0*0125 

- 0*056 

84 

- 6*4 

10.0 

78.01 

73,0 

907.7 

- 0*0097 

- 0.117 

8 5 

- 2*4 

10.0 

78,31 

70,0 

902*1 

0*4352 


- 0*129 

86 

-0 ♦ 4 

10.0 

78*74 

70.0 

903.7 

0.6694 


- 0. 132 

87 

1*6 

10.0 

80*30 

69*0 

915.4 

0.9075 


- 0.142 

88 

3.6 

10.0 

78.15 

69.0 

906,1 

1 ,1058 


- 0 . 142 

89 

5*6 

10.0 

78*51 

67.0 

906.8 

1*2202 


- 0 . 133 

90 

7*6 

10.0 

79*29 

68.0 

914.0 

1*3540 


-0 * 136 

91 

9.6 

10.0 

77*68 

67.0 

902.1 

1*5255 


- 0.153 

92 

11*6 

10*0 

80.53 

67,0 

924.2 

1.7296 


- 0.178 

93 

12*6 

10.0 

79.26 

65*0 

917.8 

1,7155 


- 0*178 

94 

13*6 

10*0 

78.44 

66.0 

910.8 

1.8233 


- 0.196 

95 

-6 • 4 

20.0 

81.20 

64.0 

930.0 

0,4258 


- 0.219 

96 

- 2.4 

20.0 

79.75 

65.0 

922.5 

0*8938 


- 0.231 

98 

9*6 

20*0 

78*67 

68.0 

914.1 

1*8714 


- 0 ♦ 255 

99 

3*6 

20.0 

77,82 

66.0 

909.8 

1.4635 


- 0.231 

100 

- 10*4 

30.0 

78*35 

64,0 

917*7 

0,3084 


- 0*288 

101 

- 6*4 

30.0 

78.57 

66.0 

913.7 

0*7973 


- 0 . 304 

102 

- 2*4 

30*0 

79,80 

66,0 

921,2 

1*2287 


- 0 . 303 

103 

1*6 

30,0 

78.09 

65.0 

907.6 

1.7413 


- 0*317 

104 

5*6 

30.0 

78.86 

64.0 

910.5 

1.9171 


- 0.311 

105 

9*6 

30.0 

78*96 

66.0 

913*2 

2,1071 


- 0.328 

106 

10*6 

30,0 

77.39 

65.0 

905*2 

2.2760 


- 0 • 349 

107 

11*6 

30.0 

79*75 

65.0 

921,8 

2.1741 


- 0 , 346 

135 

13*6 

10.0 

78*26 

79.0 

708.8 

1*7906 


- 0.125 
- 0. 134 

136 

14*6 

10.0 

79*20 

79*0 

717,7 

1.8024 


137 

11*6 

20.0 

79*00 

67.0 

716.8 

1*9553 


- 0 , 209 

138 

12*6 

20.0 

80*18 

67.0 

726.3 

1*9723 


- 0.208 

- 0.213 

139 

13*6 

20.0 

79*77 

73*0 

728*2 

1.9275 
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RUN * 
RIME 3 

ADA 

ROUGH 

FLAP 

DEF 

V ( KT ) 

T °( F ) 

76 

- 2.4 

0.0 

101.48 

73.0 

77 

- 0.4 

0.0 

102.49 

74.0 

70 

1.6 

0.0 

101.67 

73.0 

79 

3.6 

0.0 

101.06 

73.0 

80 

5.6 

0.0 

102.67 

74.0 

159 

- 10.4 

30.0 

78.76 

73.0 

160 

- 6.4 

30.0 

79.60 

73.0 

161 

- 2.4 

30.0 

78.40 

72.0 

162 

1.6 

30.0 

79.96 

72.0 

163 

5.6 

30.0 

79.28 

72.0 

164 

7.6 

30.0 

79.32 

72.0 

165 

7.6 

20.0 

79.39 

72.0 

166 

9.6 

20.0 

78.57 

71.0 

167 

8.6 

20.0 

79.10 

72.0 

160 

6.6 

20.0 

77.92 

72.0 

169 

- 6.4 

10.0 

78.30 

66.0 

170 

- 2.4 

10.0 

79.19 

66.0 

171 

1.6 

10.0 

78.61 

66.0 

172 

5.6 

10,0 

79.15 

64.0 

173 

7.6 

10.0 

79.29 

64.0 

174 

9.6 

10.0 

79.00 

64.0 

175 

10.6 

10,0 

78.22 

62.0 

176 

11.6 

10.0 

78.91 

63.0 


GLAZE 

71 

3 ROUGH 
- 2.4 

0.0 

102.00 

66.0 

72 

- 0.4 

0.0 

101.95 

73.0 

73 

1.6 

0.0 

101.88 

73.0 

74 

3.6 

0.0 

101.46 

72.0 

75 

5.6 

0.0 

100.09 

71.0 

140 

- 6.4 

10.0 

78.01 

64.0 

141 

- 2.4 

10.0 

79.61 

75.0 

142 

1 .6 

10.0 

78.88 

79.0 

143 

5.6 

10.0 

78.86 

81.0 

144 

7.6 

10.0 

78.33 

81.0 

145 

9.6 

10.0 

79.18 

80.0 

146 

10,6 

10.0 

78.24 

79.0 

147 

11.6 

10,0 

80.30 

76.0 

148 

7.6 

20.0 

79.96 

76.0 

149 

9.6 

20,0 

79.58 

74.0 

150 

8.6 

20.0 

77.88 
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